Test of constancy of speed of light with rotating cryogenic optical resonators 
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A test of Lorentz invariance for electromagnetic waves was performed by comparing the resonance 
frequencies of two optical resonators as a function of orientation in space. In terms of the Robertson- 
Mansouri-Sexl theory, we obtain (3 — S — 1/2 = (+0.5 ± 3 ± 0.7) ■ 10 -10 , a ten-fold improvement 
compared to the previous best results. We also set a first upper limit for a so far unknown parameter 
of the Standard Model Extension test theory, \{k e -) zz \ < 2 ■ 10" 14 . 

PACS numbers: 03.30.+P, 12.60.-i, 07.60.-j 
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I. INTRODUCTION 

The principle of Local Lorentz Invariance (LLI), stat- 
ing the independence of physical laws from the state of 
motion of inertial laboratories, is one of the most fun- 
damental facts about our physical world. The first pre- 
cision measurement, demonstrating LLI for light propa- 
gation, was performed by Michelson and Morley in 1887 
and its result was an essential experimental foundation 
for the advent of Relativity. LLI is incorporated as a 
fundamental symmetry into the accepted theories of the 
fundamental forces, General Relativity PJ and the Stan- 
dard Model. Numerous experiments have tested LLI 
with respect to matter and to the electromagnetic field 
and have upheld its validity until to date 0. For elec- 
tromagnetic waves the isotropy of space has so far been 
verified a the level of a few parts in 10 15 0, El H IE EI • 

New generations of tests of LLI and of other funda- 
mental symmetries (Weak Equivalence Principle, Local 
Position Invariance, CPT Symmetry) are seen as one 
important approach in the quest for a deeper under- 
standing of the forces of nature 0] . They might provide 
useful input for the development of a theory able to de- 
scribe gravity at the quantum level. In these theories, 
violations of fundamental symmetries are being consid- 
ered. Thus, the theoretical models call for improved 
experiments to either validate LLI at much higher pre- 
cision levels, or to uncover its limits of validity. 

Violations of LLI can be interpreted using so-called 
test theories. A kinematic test theory commonly ap- 
plied is that by Robertson, Mansouri and Sexl (RMS) 
Q. Here, light propagation is described relative to a 
preferred frame ("ether frame") £ in which there is no 
preferred direction and thus the speed of light cq is con- 
stant. Usually the frame in which the cosmic microwave 
background is isotropic is assumed to be this frame. 
Lorentz transformations between a laboratory frame S 
and E are replaced by general linear transformations 
which depend on the velocity v of the lab frame with 
respect to £ and on three phenomenological parameters 
a, (3, and S. These reduce to a — — ^, = |, 8 = if 
LLI is valid. In the moving frame S, the speed of light 
can be expressed to lowest order in \v\/cq as c(</>, v)/cq = 

1 + (v 2 / c o) {P ~ a - 1 + (| - P + 5) sin 2 4>) , where is 
the angle of the direction of light propagation relative 
to the velocity v of the laboratory. Thus, violations of 
the constancy of the speed of light are described by two 
nonzero parameter combinations. 

Recently, a comprehensive dynamical test theory of 



LLI violation has been developed, the Standard Model 
Extension (SME) [jj. It is based on the Lagrangian 
of the Standard Model, extended by terms that violate 
LLI and CPT. These terms describe possible violations 
in the behavior of both matter and fields, and contain a 
large number of unknown parameters, whose values can 
in principle be determined by an appropriately large set 
of (partially independent) experiments. In particular, 
the extended Lagrangian of the electromagnetic field is 
given by C = -\F^ V F^ - \(k F )^ u F^F^ . The di- 
mensionless tensor (/cf)^/^ describes violations of LLI 
and has 19 independent coefficients. Its values are de- 
pendent on the frame of reference; a frame in which 
the Sun is stationary is chosen for practical reasons. Of 
these coefficients, 10 describe polarization-dependent ef- 
fects. These can be restricted to the level of 10 ~ 32 using 
astronomical observations of the polarization of distant 
light sources [TT1 ]. 



Of the remaining 9 coefficients, one (Utr) describes 
an asymmetry of the one-way speed of electromagnetic 
waves PJ , while the others describe different aspects of 
violations of constancy, i.e. a dependence on propaga- 
tion direction and on the speed of the laboratory frame 
of reference. The corresponding eight coefficients can 
be arranged in two traceless 3x3 matrices: the antisym- 
metric (k +) ij that describes violation of boost invari- 
ance and therefore enters observable quantities weighted 
by the ratio /3® ~ 10~ 4 of Earth's orbital velocity and 
the speed of light, and the symmetric (K e _) y quantify- 
ing anisotropy of c. 



These eight coefficients can in principle be determined 
by measuring the time dependence of the resonance fre- 
quency of an electromagnetic resonator (assuming that 
particles satisfy LLI). If the electromagnetic resonator 
is stationary, seven coefficients can be determined by 
taking advantage of the rotation and orbital motion of 
Earth. Three experiments, using ultra-stable cryogenic 
resonators for microwaves and optical waves, have been 
performed recently along this line. Lipa et al. .5] took 
data for ~100 days and could constrain four coefficients 
of R e - and four linear combinations of three coefficients 
of k 0+ . Miiller et al. Q performed an experiment where 
the measurement duration was extended to a sufficient 
duration (over 1 year) that the measurement of the 7 
coefficients was achieved. Wolf et al. pfl extended the 
measurement time further and improved significantly on 
the limits. 
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The eighth coefficient, (k e -) zz , has so far not been 
determined. 
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II. METHOD 

The resonator frequency measurement involved in a 
test of LLI requires a frequency reference with a fre- 
quency instability comparable to the desired measure- 
ment accuracy 0, Q . An alternative approach is the use 
of two cavities oriented at right angles 0, . One then 
measures the difference between the two cavity frequen- 
cies, v\ (t) — Vi (t). For a two-cavity configuration, the 
beat frequency modulation can be expressed as 



2B(t)sm20(t) + 2C(t)cos2O(t), (1) 



where v\ « v 2 ~ v (2.8 • 10 14 Hz) is the average fre- 
quency and 9(t) it the angle between one cavity's axis 
relative to the south direction. Each amplitude 2B(t) 
and 2C(t) is a linear combination of eight SME coef- 
ficients weighted by time- harmonic factors. The am- 
plitude B{t) contains frequency components at 0, o>®, 
2w0, cj® ± f2® and 2u>^ ± fi®, while C(t) contains in 
addition one component at the frequency f2 ffi . Here 
is Earth's sidereal angular frequency and is Earth's 
orbital frequency. The determination of the individual 
k 0+ coefficients requires the ability to resolve the contri- 
bution of Earth's orbital motion in order to discriminate 
between frequency coefficients differing by f2®. Thus a 
measurement extending over of at least 1 year is neces- 
sary. 

The contribution of (R e -) zz to the beat frequency 
signal is 



8{ Vl {t) - v 2 (t)) 



(k e _) zz sm^ X cos29{t) + (2) 



where x is the colatitude of the laboratory. 

It is evident from eq.@ that {h e _) zz coefficient is 
only measurable by active rotation. Another advantage 
is that the apparatus has to be stable only over indi- 
vidual rotation periods, in contrast to stationary exper- 
iments, which must be stable over 12 h. Slow (com- 
pared to the rotation period) systematic variations of 
the beat frequency can be eliminated in the data analy- 
sis by fitting individual rotation periods. Furthermore, 
in a non-rotating experiment four data values can be 
obtained from a measurement lasting 12 h. 

In contrast, a rotating experiment offers the possibil- 
ity to determine two data points B(ti) and C(U) for 
every rotation (ij is the mid-time of the rotation pe- 
riod, and the period is chosen much smaller than 12 h). 
Thus, a rotating experiment offers a significant increase 
in data acquisition rate and thus a reduction of statis- 
tical noise. An overall measurement time of 1 year or 
longer is still necessary for a precise determination of 
individual coefficients. 



III. EXPERIMENTAL APPARATUS 

Crystalline sapphire optical cavities operated at cryo- 
genic temperature were used. These cavities have pre- 
viously been used to test LLI with a non-rotating setup 
as well as Local Position Invariance |T3|- Fig-D 
shows a schematic of the experiment. A mechanical 
rotation stage carries the cryostat (containing the cav- 
ities), an associated optical breadboard and most com- 
ponents (lasers, power supplies, servo systems). The 
rotation stage itself rests on an optical table that was 
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FIG. 1: The experimental set-up. Two Nd:YAG lasers are 
frequency-locked to two sapphire optical cavities located in 
a cryostat. The beams are fed to the resonators via optical 
fibers. The cavities' modes can be observed by means of two 
CCD cameras. Acousto-optic modulators (AOM) stabilize 
the power of the beams fed to the resonators. All compo- 
nents shown are mounted on a rotating table. 



not floated. The cavities were cooled to 3.4 K by means 
of a pulse-tube cooler (PTC) in which He gas is 
periodically compressed and expanded. The use of a 
mechanical cooler has the advantage of continuous op- 
eration without the need to periodically refill cryogens. 
Furthermore, slow mechanical deformations occurring 
in conventional cryostats as a consequence of change in 
cryogen levels are avoided. 

Inside the cryostat vacuum, the two cavities together 
with various optical components and photodiodes were 
mounted on a copper base plate which was supported 
by rods extending from the top flange of the cryostat. 
Copper braids provided thermal contact to the cold head 
of the cooler, and offered some compensation from the 
periodic length modulation occurring in the cooler due 
to pressure oscillation in the tubes 

The angular range of the stage was limited to about 
90 degrees by the He pressure lines connecting the cooler 
to the stationary compressor, but this was sufficient for 
the purpose of the experiment. A rotation period lOmin 
was chosen for a rotation over [0°; 90°; 0°]. Faster rota- 
tion led to disturbances. 

The cavities are high-reflection coated for wavelengths 
around 1 /im. In order to avoid rotation-induced varia- 
tions in the alignment of the laser beams with respect 
to the cavities that would result in lock errors, the laser 
waves are brought to the cavities through two optical 
fibers. Two diode-pumped monolithic Nd:YAG lasers 
were locked to the two cavities by means of a Drever- 
Hall reflection locking scheme The respective pho- 
todiodes were located in the cryostat, next to the cav- 
ities. The lasers were frequency-modulated using the 
piezo actuators glued to the laser crystals as modula- 
tors, adding the modulation frequencies to the control 
feedback signals Parts of the laser beams were 

superimposed on a fast photodiode producing a hetero- 
dyne signal at the beat frequency vi — v 2 between the 
two lasers (of the order of 700 MHz). That was mea- 
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IV. RESULTS 

The data analyzed in this publication was collected 
2 during very stable operation from 4 February 2005 to 8 
February 2005 for 76 hours. Frequency sampling was at 
6 1 s intervals. The measured tilts were time-averaged, 
weighted with the measured beat frequency sensitiv- 
~ 10 ities and subtracted from the beat data. Then, the 
~ 12 beat modulation at 1.1 Hz, caused by the PTC, was 
removed. After removal of clearly disturbed rotations, 
the remaining 432 rotation periods 9 = [0°; 90°; 0°] (la- 
belled by i) were each least-squares htted with functions 
10 a,t + 2B(U) sin26»(i) + 2C(U) cos29(t), where the coef- 
8 ficients quantifiy a (slowly varying) linear drift. 

The obtained amplitude sets {2B(ti)}, {2C(U)} are 
shown in Fig.[5] The data were analysed in the following 
ways. 

(i) Within the RMS-framework, for two orthogonal 
cavities located on Earth and considering the preferred 
4 frame mentioned above, the beat frequency can be writ- 
6 ten as in Eq. (JIJ with 

2C(t) = (1/2 -0 + S)(v 2 /c 2 )( 70 + 7i cos We r©+ 



FIG. 2: Measured 2Bv and 2CV amplitudes of spatial 
anisotropy for individual rotations, and corresponding his- 
tograms. The fit error for each data point is less than 1 Hz. 
Full line: SME model plus systematic effect. Dashed line: 
RMS model plus systematic effect. Mean values: {2B)= 2.8 
Hz, (2C)= -3.3 Hz. 



72 cos 2w e T ffi + cj i sin w e T e + a 2 sin 2w e T e ), 



2B(t) = (l/2-0 + 6)(v 2 /c 2 o )(~ /3 cosu e T (B - 



(3) 



sured with a frequency counter. 

The periodic pressure modulation (at 1.1 Hz) led to 
a residual modulation of the beat frequency between 
the two cavity-locked lasers, which amounted to ap- 
prox. 200 Hz amplitude. However, this modulation was 
later removed using Fourier filtering. The typical Allan 
standard deviation of the beat frequency (for nonrotat- 
ing resonators) was 2.5 Hz for an integration time of 
30 s, and 5 Hz (~ 2 ■ 10~ 14 ) for an integration time of 
300s, (half of the rotation period). 

The dependence of the cavity resonance frequen- 
cies on the beam powers was measured to be ~ 
8 Hz per /iW impinging on the cavities. The laser 
powers were actively stabilized using acousto-optical- 
modulators (AOM), also serving as optical isolators. At 
a typical working power of 50-100 fiW onto the cavities, 
the influence of residual power changes was thereby re- 
duced to the level of 0.1 Hz (0.4 • 10~ 15 ) for integration 
times of hundreds of seconds. 

The dependence of the beat frequency on the temper- 
ature of the baseplate supporting the cavities was on the 
order of 1.5 Hz/mK. 

Active stabilization resulted in an instability of 45 fiK 
at 300 s integration time, and thus a temperature- 
induced beat frequency instability on the order of 
0.07 Hz (2 • 10~ 16 ). 

The sensitivity of the beat frequency to tilt of the 
table supporting the apparatus was 0.06Hz//xrad (2 • 
10 _16 //^rad). The tilts along two horizontal axes were 
monitored continuously during the rotations for later 
decorrelation. Typical tilt modulation during a rota- 
tion was 50 /irad peak-to-peak, corresponding to a in- 
ferred beat modulation of several Hz. However, for the 
data analysis only the beat frequency modulation com- 
ponents at 29 are relevant; for these the inferred peak- 
to-peak amplitudes were typically less than 1 Hz. 



74 cos 2w e T e -I- cr 3 sin w e T e + cr 4 sin 2w m T m ) , 

where 70 = sin 2 x(3cos 2 — l)/4, 71 = 
(1/2) cos $ sin 20 sin 2%, 72 = cos 2$ cos 2 6 (cos 2\ - 
3)/4, 73 = cr 3 tan<I>, 74 = —0-4 tan 2$, <j\ = 71 tan$, 
02 = 72 tan 2$, cr 3 = cos $ sin x sin 29, 
04 = cos 2 cos x cos 2$ . 

T® is the time since the beginning of the data plus 
an offset that accounts for a time difference since the 
coincidence of the lab's y axis with the Y axis of the 
Sun-centered system 0. The direction of the Sun's ve- 
locity v relative to the cosmic microwave background 
is given by the right ascension $ = 168° and the dec- 
lination = —6°. In fitting the above expression to 
our data, we take into account that the finite average 
values of 2B and 2C are likely to be due to a (con- 
stant) systematic error. We model this by adding a 
contribution b sys t sin 2(9 + c syst cos 29 to Eq.Q), and fit 
b syst , c S y S t > — S — |. Effectively, this means that only 
the modulation of the {2B, 2C} amplitudes by Earth's 
rotation contributes to the fit result for — S — \ . 

The fit yields - S - \ = (+0.5 ± 3) • 10" 10 and 
is shown in Fig.|3 Due to the experimental error in 
the determination of the tilt sensitivities there is an ad- 
ditional (systematic) error of ±0.7 • 10~ 10 . Our result 



TABLE I: Fourier amplitudes determined from the experi- 
ment. All quantities are in units of 10 -16 . 





Basis 


Value 


Statistical error 


Systematic error 


Co 


1 


-59 


3.4 


3 


Ci 


sin (cj©T©) 


-3 


1.5 


0.5 


c 2 


cos (w©T e ) 


11 


2 


0.5 


C 3 


sin(2w©r e ) 


1 


2 


0.5 


c 4 


cos(2a) e T e ) 


0.1 


2 


0.2 
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is about factor 10 lower than the previous best results 
(-1.2 ± 1.9 ± 1.2) • 10" 9 19] and (-2.2 ± 1.5) • 1CT 9 0. 

(ii) Using all our data we can obtain a fit of the 5 
Fourier amplitudes of 2C(t*) at the time t* ~ 6 February 
2005, see Table H] These amplitudes are linear combina- 
tions of the and k 0+ coefficients, where the weights 
of the latter depend on Earth's orbital phase (see Ap- 
pendix E of 9J). For the fit we make use of the approx- 
imate relationship between the functions B(t) and C{t) 

20] and we include a coefficient b in B(t) in order to de- 
scribe a systematic effect. Note that due to the limited 
extent of the data over time, the frequencies + fi^g 
and W0 — J7m cannot be distinguished and therefore it 
is not possible to extract the individual coefficients k e - 
and k 0+ from the Cj coefficients. 

(iii) The results of the microwave cryogenic exper- 
iment found the elements of (n e -) (except for 
(k e S) zz ) and the elements of (3q (k +) to be zero within 
several parts in 10 -15 . If we assume these elements to 
be zero, we can use only the Co-coefficient to determine 
{h e -) zz , i.e. Eq.|(5J| is truncated at the first term. The 
result is {k e ^) zz = (-2.0±0.2) • 10~ 14 . As this nonzero 
value is likely due to a systematic effect, we may state 
that the value of \(k e -) zz \ is probably less than 2-10~ 14 . 

We believe the nonzero averages of the 2B and 2C am- 
plitudes are due to a systematic effect of thermal origin. 
A (small) thermal gradient is present in the laboratory, 
through which the apparatus rotates. The magnitude of 
the effect is consistent with our measured temperature 
sensitivity and the temperature modulation amplitude 



measured at the top of the apparatus. 

Overall, our experiment is limited firstly by the slow 
maximum rotation speed allowed by the system. A 
faster, variable rotation speed would in principle allow a 
more detailed study of the systematic effects or a partial 
suppression thereof. However, faster rotation could not 
be used with our setup. Secondly, the laser frequency 
lock stability is limited by the small signal-to-noise ra- 
tio (low cavity throughput), which limits the ability to 
more precisely characterize the systematic effects. 

In conclusion, we have described a Michelson-Morley- 
type experiment that has led to the first determination 
of a limit for (k e ^) zz of the Standard Model Extension, 
\k zz \ < 2 • 10~ 14 . An analysis of the data in the RMS 
framework yielded an upper limit for a hypothetical di- 
rection dependence of c an order of magnitude lower 
than the best previous measurements. We note that 
the present results were achieved using a measurement 
duration one hundred times shorter. The accuracy of 
the experiment appears to be limited by a thermal sys- 
tematic effect. A significant improvement of our results 
would require major modifications to the apparatus. 
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